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Abstract. We introduce bi-variate and multivariate force of mortality functions. The pattern of mor- 
tality in a population is one of the strong influencing factors in determining the life expectancies at 
^ ' various ages in the population. The reasons behind declining forces of mortality could be studied using 

the proposed functions. 
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1. Introduction 

^ One variable force of mortality, with respect to age, x of an individual is one of the central 

■ topics of study in the actuarial mathematics |SmithAMM1948| . It is often termed as instantaneous rate 
of death at an age x. Suppose, forces of mortality is measured on two variables {x,y), (one being age, 
other can be some influencing factor on mortality), then if we plot the force of mortality fi{x,y) on the 
xy— plane, then for an arbitrary point (xo,?/o); we can write. 



;i.l) ^^{xo,yo) = / / fi{x,y)dp 




n 



where is a region such that {xq, y^) G and p(fi) is area of the region. Since yU is continuous in the 
univariate case, if we assume the same holds for /x (x, y) in f2, it has an upper bound say, /ii and a lower 
bound say, /iq in the region such that 
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Multivariate force of mortality functions can help in better understanding of future longevity and 
causes of decline in mortality rates. There are studies which consider mortality decline or longevity pro- 
jections of humans with respect to age only (for example, see |TuljaLiBoeNature2000 , IWilmothSciencel998] ). 



Such studies can be handled using univariate analysis of standard force of mortality functions. However 
such studies can be extended to incorporate several variables that can explain decline in mortality rates 
or causes of increase in longevity using multivariate forces of mortality functions. Further analytical 
properties such as rate of change in forces of mortalities of bi-variate force of mortality functions are 
found in the next sections. We will begin with two examples and we illustrate their numerical properties. 

Example 1. Suppose l{x, y) = 1 — x'^y''^^ /K for a,b,K E Z+. See the definitions of fix{x, y) 
and fi{x,y) in eq. (13. 6p and eq. (13.71) in the section 3. We have, 

, , ax^'-Y^'^ 
^^x{x,y) -- 



fiy{x,y) 



^^ay{b/VK)-l 



K{K - x^y^/^) 



Example 2. Suppose l{x,y) = a^b^'^ for 10 < x < 80, 1 < ?/ < 5, a = 5, 6 = 7. jj,x{x,y) = — 3x^log(6) 
and y) = - \og{a)/2y/y. 

Force of mortality functions also associated with continuous life table functions |SmithAMM1948] . 
Life table is a mathematical model describing how individuals born at same time (for example, a cohort 
of new born babies) survive over the years at various ages until the last individual dies. Life table is 
constructed based on present or past mortality pattern (i.e. mortality rates at each age, (say x) per 
fixed number of individuals in the same age x per year or for the year (0,t)) in the population and 
assumed that this pattern will remain the same until the individual at last age dies. This table can 
be used to construct synthetic population at each age x at time t or for time interval (0,t), it do not 
have mechanism to take care of future changes in the mortality pattern after t. Suppose l{x) denote the 
number of individuals at age x out of l{0) newly born individuals in a life table with continuous partial 
derivatives up to order (A; + 1). Then the number of individuals at ages, x + Ax and x — Ax are denoted 
by l{x + Ax) and l{x — Ax) can be obtained from Taylor series expansion evaluating at xq as 



/(xo + Ax) = /(xo) + Ax/'(xo) + (Ax)2^p^ + (Ax)^^l^ + ... 



,j('=)(xo) , r (Ax)'^/('=+i)(t) 



(1.2) +<^-)^'^ + i / 
l{xa - Ax) = l(x„) - Arf(x„) + (Axf'-pl - (Axf'-pl + ... 

(1.3) +(_l)^(A.)'?^+r(M!^* 

Here x is the fixed age between and u, the maximum age of life. The number of survivors at age 
X with some other relevant factor (for example, marital status, education level, climate, food habits. 
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Figure 1.1. l{x,y) in Example 1. For various combinations of values of a and b we 
have drawn l(a — 1) to l{i — 1) by fixing K = 100, 10 < s < 80 and 1 < y < 5. 
Following are the combinations of a and b for each figure: l(a — 1) : a = 0.1,6 = 0.9, 
1(6 - 1) : a = 0.2,6 = 0.8, l(c - 1) : a = 0.3,6 = 0.7, l(ci - 1) : a = 0.4,6 = 0.6, 
l(e - 1) : a = 0.5,6 = 0.5, 1(/ - 1) : a = 0.6,6 = 0.4, 1(^ - 1) : a = 0.7,6 = 0.3, 
l(/i - 1) : a = 0.8, 6 = 0.2, l(i - 1) : a = 0.9, 6 = 0.1. 
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l(a-2) l(d-3} l(g-2) 





Figure 1.2. fix{x,y) in Example 1. For various combinations of values of a and b we 
have drawn l(a — 1) to l(i — 1) by fixing K = 100, 10 < a; < 80 and 1 < y < 5. 
Following are the combinations of a and b for each figure: l(a — 1) : a = 0.1,6 = 0.9, 
1(6 - 1) : a = 0.2,6 = 0.8, l(c - 1) : a = 0.3,6 = 0.7, l(ci - 1) : a = 0.4,6 = 0.6, 
l(e - 1) : a = 0.5,6 = 0.5, 1(/ - 1) : a = 0.6,6 = 0.4, 1(^ - 1) : a = 0.7,6 = 0.3, 
l{h-l):a = 0.8, 6 = 0.2, l(z - 1) : a = 0.9, 6 = 0.1. 
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l(a-3) l(d-3} l(g-3) 




Figure 1.3. fiy{x,y) in Example 1. For various combinations of values of a and b we 
have drawn l(a — 1) to l{i — 1) by fixing K = 100, 10 < a; < 80 and 1 < y < 5. 
Following are the combinations of a and b for each figure: l(a — 1) : a = 0.1,6 = 0.9, 
1(6 - 1) : a = 0.2,6 = 0.8, l(c - 1) : a = 0.3,6 = 0.7, l(ci - 1) : a = 0.4,6 = 0.6, 
l(e - 1) : a = 0.5,6 = 0.5, 1(/ - 1) : a = 0.6,6 = 0.4, 1(^ - 1) : a = 0.7,6 = 0.3, 
l{h-l):a = 0.8, 6 = 0.2, l(z - 1) : a = 0.9, 6 = 0.1. 
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geographic region etc) y evaluated at (a;o,l/o) can be obtained form two variable Taylor expansion. 
Assuming continuous partial derivatives for / [xq + Ax, Uq + Ay) up to order 3, we can the expansion of 
two variable survival functions as follows: 

{ dl 1 { dl 1 (Ax)2 f 
l{xo + Ax,yo + Ay) = / (xq, ?/o) + Ax | — (xq, yo) | + Ay | — (xq, yo) | H 2!~ | Sx^ 

[xo, yo) > + AxAy <^ — — (xo, yo) > + < tt^ [xq, yo) 



2! [dy"^ J [dxdy J 3! [dx^ 



3! {dy^ ' j 2 [dx'^dy 

,1.4) +^J^|_(,„,,„)| + + 

[9/ 1 f 9/ 1 (Ax)^ f 9^/ 
Z (xo - Ax, 2/0 - Ay) = / (xo,?/o) - Ax | — (xo,?/o)| - Ay |— (xo,yo)| H 2!~ | 9x2 

(Ay)2r52/, ^ ^ { on , A (Ax)3r93/, 

(a;o, l/o) \ + AxAy <^ — — (xq, yo) \ ^ { [Xq, Vo) 



2! [ 9y2 J 1^ 5x9y J 3! [ c?x^ 

\2 



(Ay)3 ( dH ^ ,1 (Ax)' Ay f S^/ ^ ^] 
-^|5^(^o,l/o)j 

\ dxdy"^ 



2. Analysis of first order equations 

We try to analyze univariate force of mortality functions by considering higher order derivatives are 
continuous. Assuming x — )■ Xo, in eq. dOD and eq. ([LS]), we will have f!^ {{Ax)H^^+^\t) / k\) dt 0. 
Suppose f'^"'\xo + Ax) and /*^")(xo — Ax) denote Taylor expansion equations when ignoring the (ri + 1)*'* 
derivatives and beyond for n = 2, 3, ... in eq. (II. 2p and eq. (11.31) . then by sequentially ignoring the terms 
beginning from the term ^^^^ /^"^(xo) in eq. (II. 2p and eq. (11.31) for n = 2, 3, ... , we will obtain following 
equations: 

(2.1) /(i)(xo + Ax)-/(^)(a;o-Ax) = 2Axr(xo) 

(2.2) /(3)(xo + Ax)-/(3)(xo-Ax) = 2Axl' (xo) + '^^^^l^'\xo) 
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f^'\xo + Ax)-f'\xo-Ax) = 2Axnxo) + ^^l^'\xo) + ^^l^'\xo) 



3! 



5! 



/(^'=-^)(xo + Ax) - f<''-'\xo - Ax) = 2Axr(xo) + ^-^f'\xo) + ... + ^-^^^^^^"'''^^o) 



(2.3) 



Therefore, 

E •=! [{f^'^'K^o + Ax) - /(2^+i)(xo -Ax)}- {f'^-'Kxo + Ax) - f^='^\x, - Ax)]] 



(2.4) 



2(Ax) 
3! 



'>\xo) + ^-^l^'\xo) + ... + ^-W^l'''^'\^o) 



5! 



{2k + l)\ 



Hence we will obtain 



E-=i [{/^'''^'Ha^o + Ax) - /(2^+i)(xo - Ax)} - {/(2^--i)(xo + Ax) - /(2^-i)(xo - Ax)}] 



(2.5) {f^''^'\xo + Ax) - f^'^'-'\x, -Ax)]- {/«(xo + Ax) - /«(xo - Ax)} 

Using the relation I'^'^^xq) = [f^^\xo + Ax) — /'■^•'(xq — Ax)] , we write, 



dl{x) 



dx 



x=xo 



2Ax 

k 



|/(2fc+i)(a;„ + Ax) - f'^^''+'\xo - Ax)] 



(2.6) Yl K/^''^'H^o + Ax) - f^''+'\xo -Ax)]- {f^'^~'\xo + Ax) - f'^''-'\xo - Ax)]] 

i=i 

Dividing the eq. (12. ip by the term 2Ax on both the sides and integrating it from age x to x + m, we 
obtain. 



x+m 



{fW{y, + Ay)-fW{yo-Ay)] 
2Ay 



dy 



x+m 7 

dy 



dy 



y=yo 

x+m ^ ^ 



l{y) dy 



y=yo 



(2.7) 



x+m 



M\y=yoKy)\y=yody 



where, yu(x), the force of mortality function, which is defined as (— l//(x)) {d/dx)l{x). 
Dividing the eq. (12. 2p by the term 2Ax on both the sides and integrating it from age x to x + m, we 
obtain. 



(2.8) 



r-x+m 



{/(3)(yo + Ay)-/W(yo-Ay)} (AyY^ 



2Ay 



3! 



i^'\yo) 



dy 



x+m 



Ky)\y=y, Ky)\y=y,dy 
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Now, multiplying and (— to the eq. fl2.6p . and integrating between ages x and 

X + m, we will obtain, 



x+m 



1 d 

l{y) dy 



Ky) 



y=yo 



x+m 



2Ay 



{f^''^'\yo + Ay)-f^''^'\yo-Ay)} 



{f'^'^-'\yo + Ay)-f^'^-'\yo-Ay)}]] 



[2.9) 



x+m 



dy 



Ky) 



dy 



y=yo 



ml'iyo 



If there are any deaths during the age {xq + Ax) to (a;o — Ax), then f^^^ 
In the absence of deaths, we have /^^-"(xq + Ax) = f^^^^XQ — Ax) and 
true for difference of other higher order expressions. 



(xo + Ax) </«(xo-A 



X 



< 0. This argument is 



x=xo 



3. Analysis of second order equations 

Suppose g^'^\xo + Ax, yo + Ax) and g^"'\xo — Ax, yo — Ay) denote equations when ignoring the terms 
from {n + 1)*'^ partial derivatives and beyond for n = 2,3,... in eq. (11. 4p and eq. (11. 5p . then by 
sequentially ignoring the terms beginning from the n^^ order partial derivative terms in eq. (11.41) and 
eq. (ll.Sp for n = 2, 3, ... , we will obtain following equations: 

g^'^ (xo + Ax, yo + Ay) - g^''^ (xq - Ax, yo - Ay) = 2Ax 1 1^ (xq, yo) | + 2Ay | ^ (xq, yo) | 
g^^^ (xo + Ax, yo + Ay) - g^^^ (xq - Ax, yo - Ay) = 2Ax <^ — (xq, yo) ^ + 2Ay <^ — (xq, yo) 



dx J {dy 

(Ax)3 ( OH , A (Ay)^ ( OH , 



{Ax)^Ay ( dH ^ 



2 

\2 



Ax (Ay)' ( d'H ^ 
2^19^9^^^°'^°) 
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(a;o + Ax, yo + Ay) - g^'''-'^ [xq - Ax, y^ - Ay) 



(2fc-l) 



2Ax { ^ (xo, yo) I + 2A|/ 1 ^ {xq, yo) | + 



2 {Ax 



,2k-j-l 



(Ayy 



{2k-l)\ 



2k-l 
.j=0 



2k - 1 



{xo,yo) 



Therefore, 



E [{9^'"^'^ (^0 + + - 9^'''''^ (^0 - Ax, yo - Ay)} 

-{9'^''~'Hxo + Ax,yo + Ay)-g'^"^~'Hxo-Ax,yo-Ay)}] = ^|S(xo,2/o) 



3! (9x3 



(Ay)'^ r 



3! 1^ dy 



{xo, yo) 



(3.1) 



+ 



2 (Ax 



,2fc-i+l / A„Ai 



{Ayy 



i2k + l)\ 



(Axf Ay ( dH , 

'2k+l 

-i=o 



Ax (Ay)^ f a^/ 



j I [ ax(2fc-J+i)52/^ 



(a;o, 2/0) 



Force of mortality for the two variables (x, y) is evaluated at the point (xq, yo) using partial derivatives 
as follows: 



(3.2) 

(3.3) 

where, we define, 
(3.4) 

(3.5) 

(3.6) 

(3.7) 



dfJ^ jx, y) 
dx 

dfJ, jx, y) 
dy 

/i(xo + Ax,yo) 
l^{.xo,yo + Ay) 
fJ'x {xo,yo) 
^■y {xo,yo) 



lim 

Ax-5>0 



lim 

Ay-i>0 



/X (xq + Ax, yo) - /i {xo^yo) 
Ax 

{xo, yo + Ay) - jj, {xo,yo) 
Ay 



1 



dl (xo + Ax, yo) 



I (xo + Ax, yo) dx 

1 dl (xq, yo + Ay) 

I {xo, yo + Ay) dy 
1 dl{xo,yo) 



I (xo, 2/0) dx 

1 dl{xo,yo) 
I (xo, 2/0) ^1/ 



for eq. ([S3) 
for eq. (IX^ 



Using eq. (11. 4p and as Ay — i- 0, we obtain. 
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l{xo + Ax,yo) = l{xo,yo) + Ax {—{xo,yo) > ^ ^ ^ (^o, Z/o) ^ H ^ ^ (a^o, l/o) 



(9x ' j 2! l^(9x^ ' j 3! |^(9x'^ 

(3.8) 

Using eq. (II. 5p and as Ax — i- 0, we obtain, 

. . i dl , A (Ay)^ ( dH , (Ayf ( dH{xo,yo) 

l{xo,yo + Ay) = U^o,yo) + Ay|-(xo,yo)| + ^|^(xo,yo)| + ^|— 

(3.9) 

Therefore, 

dl ^ A N f dl , A ^ f dH , ,] (Ax)^ ( an ^ A {Axf ( dH ^ 

-(xo + Ax,yo) = |a^(^o,yo)|+Ax|— (xo,yo)j + ^|^(xo,l/o)| + ^|^(xo,l/o) 

(3.10) 

9h A^ f 91 , A ^ f d^h A (Ay)^ f dH ^ A {Ayf f d^l ^ 

— (xo,yo + Ay) = |_(xo,yo)|+Ay|— (xo,yo)| + ^j-|^(xo,l/o)| + ^p|^(xo,Z/o) 

(3.11) 

Let us now derive the equation of the type eq. (13.11) with the conditions Ay — )■ and Ax — )■ and 
extending these equations up to the general term. Suppose, Ay — )■ in the equations (II. 4p and (II. 5p , 
then 

. . i 91 , A (Ax)2 ( d^l , 
l{xo + Ax,yo) = l{xo,yo) + Ax < — {xo,yo) > ^ ^^(^o,yo) 



dx I 2! 1 dx"^ 

. . i dl , A (Ax)2 ( dH , 
l{xo-Ax,yo) = l{xo,yo) - Ax < — {xo,yo) } ^ ^ tt^ (a^o, J/o) 



dx J 2! [ dx^ 

(3.3) „„,,„)}+... 

and Ax — )■ in the equations (II. 4|) and (II. Sp . then 

( dl ] {AyY f d'^l 
l{xo,yo + Ay) = / (xq, yo) + A?/ 1 — (xq, yo) | H 2P | (^O' ^o) 1" + 



. . ] dl , A (Ay)' i dH , 
l{xo,yo-Ay) = / (xo,yo) - A|/ <; — (xo,yo)| H 2)~ ^^^'^^^ ' 



(3.15) +(-l'"^{P<^°-''»)^ + " 
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Sequentially, ignoring the n*^^ order terms from eq. fl3.12p and eq. fl3.13p . and denoting these new 
equations as g^^ {xq + Ax, jjo) and g["'^ {xq — Ax, yo) for n = 2, 3, we obtain below set of equations. 

f dl 

(3.16) g'i^\xo + Ax,yo) - gi\xo- Ax,yo) = 2Axl — {xo,yo) 

&_ 

dx ^ j ' 3! 1^ dx'^ 

gf^~^^ (xo + Ax, yo) - g[^^'^^ (xq - Ax, yo) = 2 Ax \ ^ (xq, yo) j + ^^^j^ | ^ {xq, yo) \ + 

+ 



«f'(.. + A..,.)-.f)(.„-A..,.) ^ 2A.^#(.„.,„)UH(^m,.„,,„) 



dx ^ j 3! 1 dx^ 
{2k -1 



Therefore, 

k 

J2 [{91'^'^ (^0 + Ax, yo) - gf'^''^ (xo - Ax, yo) } - [gf'"'^ {x, + Ax, y,) - gf^''^ (xo - Ax, y,) 



j=0 



A2(AX)(2^-+1) [9(2^-+!)/ , ,1 

= js-f ''^^^ (a^o + Ax, yo) - 9?''^^^ {xo - Ax, ?/o)} - js'i^^ (a^o + Ax, yo) - fi-f ^ (xo - Ax, yo) 

Sequentially, ignoring the n*^^ order terms from eq. fl3.14p and eq. fl3.15p . and denoting these new 
equations as g^^ (xo, yo + Ay) and g^^ (xq, yo + Ay) for = 2, 3, we obtain below set of equations. 

(3.18) 5-5^^ (xo,i/o + Ay) - ^fj^^ (xo,yo - Ay) = 2Ay|— (xo,yo) 

dy J ' 3! 1^ dy^ 

9?'^^ (xo, yo + Ay) - gf^'^^ (xo, yo - Ay) = 2Ay \ ^ (xo, yo) j + ^^^f^ 1 14 (^0, Z/o) + 



(3) . . X (3) . .X . ) 5^ . 2(Ay)3 ( d^l , 
92 {xo, yo + Ay) - g'2 (xo, yo - Ay) = 2Ay <^ — (xo, yo) > H — < tt^ (xq, yo) 



dy } 3\ { dy 

2(Ay)(2'=-i) r 9(2^-1)/ 
+ (2A:-1)! 1^75^^^°'^°^ 



Therefore, 



E [{92'^'^ (^0, 2/0 + Ay) - (xo, yo + Ay) } - [gi''-'^ (xo, yo + Ay) - gi'^''^ (xo, yo + Ay) }" 



j=0 
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= (xo, z/o + Ay) - P^^^^ (xo, + Ay) } - {/(I) (xo, yo + Ay) - /(^) (xo, yo - Ay) } 

Now substituting the eqs. (I3.8p and (l3.1Qp in the eq. (13 ■4p . we get 



^ (Ax)^- r d^^H 



(3.20) u(xo + Ax,yo) = — f i y ., i r^ ■ , , 

[!(:^o.!/o) + E~.4^{ii(^<..!/»)}]^^ 

and, substituting the eqs. (13.91) and (I3.1ip in the eq. (13. 5p . we get 



(3.21) (xo, 2/0 + Ay) 



U^o,2/o) + Er=i^{0(^o,l/o) 



^ (Ay)^- r 



Using eq. (I3.16P and eq. (I3.18p . we obtain forces of mortahties for two variables as follows: 

[gi^ (xo + Ax, yo) - g^^ (xq - Ax, yo) | 



At (xo.yo) 



2Ax/ (xo, yo) 
[g? {xo, yo + Ay) - g^^'^ (xq, yo - Ay) | 



2 Ay/ (xo,yo) 

Hence the derivative forces of mortalities are as follows: 



dfi (x,y) 
dx 

(3.22) 

dfi (x,y) 



dy 



(3.23) 



lim — — 

Ax-^O Ax 



lim — — 

A{/-s>o Ay 



g^^^ (xo + Ax, yo) - gi' (xo - Ax, yo) 
2Ax/ (xo,yo) 



(1) 



oo (Axy an 
i=i j\ &J 



{xo,yo) 



\^g2^ (a;o, 2/0 + Ay) - g'i'^ (xq, yo - Ay) 


1 v-oo (Ay)i J 


[ ay^+i 


^a;o, yo) 




2Ay/ (xo,yo) 


^ (a;o,2/o) + EJli 


{^yV J 
i! 1 


(a;o,2/o) 





4. Analysis of third order, equations 

Suppose s be the function number of survivors at age x with two more influencing variables y and 
then the three variable survival function s (x, y, z) evaluated at {xq, yo, Zq} can be written as 
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s (xo + Ax, yo + Ay, zq + Az) 



(4.1) 

Similarly, 



s {xo,yo,Zo) + Ax 1^ (xo,2/o,^o)| + {xo,yo,Zo) 

( ds ] Ax^ f d^s 

+Az < — {xo, yo, zo) > + < ^ (xo, yo, zo) 

Ay^ (dh , A Az^ (d^s , 

+AxAy <J ^TTTT- (a^o, 1/0, 2:0) > + AxAz < (xq, yo, zq) 



+AyAz 



n=k 



dxdy 
dydz 

E 

.ni,n2,ra3 



(xo,yo,2o) ^ + 



1 



ni\n2\n^\ dx'^^dy^^dz 



"3 



Ax^^Ay^^Az 



s (xo - Ax, yo - Ay, zo - Az) 



(4.2) 



s (xo, I/O, 2^0) - Ax I — (xo, yo, zo)>-Ayl— (xq, yo, zo) 
rSs , ,1 Ax2 \dH , 

Ay2 rs^s , ,1 Az^ (d^s ^ 

+AxA?/ <j ^^-^ (xo, z/o, ^0) \ + AxAz <^ (xo, yo, zo) 



+AyAz 



dxdy 
dh 
dydz 



(xo, 2/0,^0) (■ + 



n=k 



.ni,n2,ra3 



1 



d''s 



nilrislng! dx'^^dy'^^dz'^^ 



Ax"^Ay"2A2' 



Here = rii + n2 + ns. Let /i*^"^ (xq + Ax, yo + Ay, zo + Az) and /i*^") (xq — Ax, yo — Ay, zo — A 
denote equations after ignoring the terms with derivatives beginning from the (n+1)*'^ order [n = 1,2,. 
in the eqs. (14.11) and (14. 2p . We will obtain following equations: 



(xo + Ax, yo + Ay, zo + Az) 
-h^^^ (xo - Ax, yo - Ay, zq - Az) 



ds 1 r ds 

2Ax { — (xo, yo, Zo) | + lAy | ^ (^0, Vo, zo) 



ds 



+2A^ { ^ (xo, yo, Zo) 



MULTIVARIATE FORCE OF MORTALITY 



1) (xo - Ax, I/O - A|/, 2;o - A^;) J t dx J J 

+2A2; (a;o,yo,^o)| + 

h V Ax^iAy''^^ 

^ ni!n2!n3! (9a;"i(9v"292;"3 

ni,n2,n3 

(here 2k — 1 = rii + n2 + n^) 



Therefore, 



^ (xo + Ax, yo + Ay, zo + Az) - h^"'+'^ {xo - Ax, y, - Ay, zo - Az)} 



(a;^ + Ax, yo + Ay, ^0 + A^) - ^(^^-i) (3;^ - Ax, yo - Ay, + A^) }] = 



2 i9^s 

Ax"^ Ay"-'^ Az^^ 

ni,n2,rt3 

(here 3 = ni + n2 + na) 



2 /)(2fc+i) „ 
= + •••+ y , , — Ax"iAy"^Az"3 

m,n2,n3 

(here 2A; + 1 = ni + 712 + ris) 
We define following three force of mortality functions evaluated at (xq, yo,Zo): 



IJ'x {xo,yo, Zo) 
l^y {xo,yo, Zo) 
l^z {xo,yo,zo) 



1 


ds (xo,yo,^o) 


s (xo,yo,^o) 


9x 


1 


9s (xo,yo,2;o) 


s (xo,yo,^o) 


dy 


1 


ds (xo,yo,^o) 


s (xo,yo,^o) 


92; 



and further we define three functions of forces of mortahty as follows: 
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1 ds{xo + Ax,yo,zo) 



fi{xo,yo + Ay,zo) 
/i (xo, 2/0, ^0 + Az) 



Q ( T* ^ 1 /\ /J 'y^ 1 

6 l^xo -f- ^x, yo, zq) 


dx 


1 


ds {xo,yo + Ay,zo) 


s (a;o,|/o + Ay,zo) 


dy 


1 


ds {xo,yo,zo + Az) 


s (xo,yo,2;o + Az) 


dz 



Using the above definitions, we obtain following rates evaluated at (a;o, yo, Zq): 



dfi {x, y, z) 


= lim 

Aa;->0 


/i (xo + Ax, yo, zo) 


- fi{xo,yo,zo) 


dx 


Ax 




dfi (x, y, z) 


= lim 

Aj/-s>0 


(xo, yo + Ay, Zo) 


- fi{xo,yo,zo) 


dy 


Ay 




dn (x, y, z) 


= lim 

Az-¥0 


n{xo,yo,zo + Az) 


- fi{xo,yo,zo) 


dz 


Az 





By taking pairs of limits(A?/ — )■ 0, Az — t- 0), (Ax — 0, Az — )■ 0), and (Ax — )■ 0, Ay — > 0), separately 
in the equation (14. ip . we obtain following three equations: 



ds ^ Ax" ( d^s 

s{xo + Ax,yo,zo) = s (xQ,yQ, zq) + Ax { — {xo,yo, zo) > -\ ^ | ^x" 



(4.3) 



ds 1 Ay" f(9"s 

s{xo,yo + Ay,zo) = s{xo,yo,zo) + Ay { —{xo,yo,Zo)\ -\ ^ |^ ^O' + 



(4.4) 



s{xo,yo,zo + Az) = s (xq, 2/0, ^o) + A^ (xq, yo, ^o)| H ^ {f^ ^°''} 

(4.5) 

By taking pairs of limits(A?/ — t- 0, Az — )■ 0), (Ax — t- 0, Az — )■ 0), and (Ax — )■ 0, Ay — t- 0), separately 
in the equation (14. 2p . we obtain following three equations: 

( ds ^ Ax" f 1 
(4.(a)o - Ax, yo, Zo) = s (xo, yo, zo) - Ax <^ — (xo, yo,zo)> ^ h ("l)"^^ \ ^ (^0, 2/o, ^0) ^ H 

f 9s 1 Ay" f9"s 1 

(4(^:^0, ?/o - Ay, Zo) = s (xo, yo, ^o) - Ay | — (xq, yo, 2^0) | + ■ ■ ■ + | ^ (^0, ?/o, 2^0) | + ■ • • 

( ds 1 A;z" f 9"s 1 
(4.#o, 2/0, ^0 - Az) = s (xo, yo, 2:0) ~ ^^ig^ (^O' 2/o, 2:0) > H h ("1)"-^ j ^ (^0, l/o, zo) > -\ 
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Suppose h'l (xo + Ax, yo, zq), (xq - Ax, yo, zo); (xq, yo + Ay, zq), (xq, yo - Ay, 2:0) and 

/ig (xq, yo, Zq + Az) , /ig (xo, yo, Zo — Az) for n = 1,2, 3, ... denote the functions by ignoring the terms 

from the order {n + 1) in the pairs of equations (??), fl4.6p ; fl4.4p . (14. 7p and ( 14. 5p . (14. 8p . then we will 

obtain following three series of sequences of difference functions: 

^[^^ (xo + Ax, yo, Zo) - h^i^ (xq - Ax, yo, Zo) 
hi^ (xo + Ax, yo, Zo) - hf^ (xq - Ax, yo, zo) 



2Ax 



ds 
dx 
ds 



(xo, yo, Zo) 



2Ax <{ — {xo,yo,zo) )■ + 



2Ax^ r d^s 



3! 



dx^ 



{xo,yo,zo) 



hf" (xo + Ax, yo, Zo) - hY" (xq - Ax, yo, Zo) 



{2k-l) 



( ds 1 2Ax^ ( d^s 

2^a; I — (xo, yo, zo) | + | ^ (xq, yo, zq) } + 

2^^(2fc-l) J ^(2fc-l) 

(2fc-l)! la^M^"^"'^"'^"^ 



Therefore, 



k 

E 



/if^'^^^ (xo + Ax,yo,zo) - hf^^^'^ (xq - Ax, 1/0,^0)} - 
/^S^^"^^ {xo + Ax, yo, Zo) - hf^'^^ (xq - Ax, yo, zo) } 



2Ax=' r d^s 



(4.9) 



/i2^^ (a;o, 2/0 + Ay, zo) - h)^'' (xq, yo - Ay, zq) = 2Ay { ^ (xq, yo, Zo) 



{xo,yo,zo) } + 
+ (2A: + 1)! 



3! 1^ dx^ 



,(1) 



9s 



,(3) 



9?/ 

ds 



2Ay^ f S^s 



3! [ dy 



(xo, 2/0,^0) 



(a^o, 2/0 + Ay, Zo) - hf'' (xo, yo ~ Ay, zo) 



,ds , ,} 2Ay3 r 

2Ay <^ — {xo,yo,zo) } + 4 



+ 



dy 
2Ay2'=-i 



3! [dy 



{xo,yo,zo) 



{2k-l)\ \dy 



d 



2fc-l , 



,2k 



- {xo,yo,zo) 



Therefore, 
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k 

E 



(xo, z/o + Ay, zo) - h^2^~^^ {xq, yo - Ay, zq)^ 



2 At/3 ( d^s 



(4.10) 



{xo, yo, Zo) } + 



3! t dy^ J 



h'i^ {xo, yo, Zo + A^;) - {xq, yo, zo - Az) = 2Az | ^ (xq, 2/0, 2:0) | 



hs\xo,yo,zo + Az) - h^^' {xo,yo,zo- Az) = 2Az <j |- (xq, yo, ^0) + 



as 



2Az^ (d^ 



3! 



(xo,yo,^o) 



h'i'' (xo, yo,zo + Az) - hf^ {xo, yo, ^0 - Az) 



{ ds "1 2Ar'^ fa^s 

2 A^ I ^ (a^o, ?/o, ^0) j + I (^0' 

2A2(2fc-l) J^(2fe-l) 

(2ifc-l)! (a^M ^^0'^°'^°^ 



Therefore, 



k 

E 



/if ^"^^^ (xo, yo, Zo + A2;) - h^^^""^^ {xo, yo, zo - Az) j 



,(2i+i) 



|/if ^ {xo, yo,zo + Az) - /if ^ (xo, yo, ^0 - A^) | 



2Az^ (d^ 



{xo,yo,zo) > + 



(4.11) 



3! t dz^ 

2A^(2fc+l) fd{2k+l)^ 



+ 



{2k + l)\ 



The rate of changes in the survival function with respect to one variable and corresponding forces of 
mortalities for three variables can be obtained using the following derivations. 



ds(xo + Ax,yo,zo) ds(xo,yo,zo) , a ^ f _ n1 , , Aa;"f(9"+^s 



dx 



dx 



Ax" f (9"+^s 1 
+ ^ ^ (^0, ^0, ^0) j> + ■ ■ • + ^ I {xo, yo, Zo) I + 



ds {xo, yo + Ay, Zo) ds (xo, yo, zo) . ( d^s \ Ay'' ( a^+^s , 
dy = d-y + ^^|a^(^o,yo,.o)j + --- + ^|^(xo,yo,.o)K 



ds {xo,yo,zo + Az) 
dz 



ds (xo yo Zo) f d^s 1 Az''^ f (9"~'~^s 

+ ^^\-Q;^i^o,yo,zo)^ + --- + -^^-g^{xo,yo,zo)} + 



dz 



H {xo + Ax, yo, Zo) 



s {xo, yo, Zo) + ET=i (0 yo, Zo)}] 



i=0 



Ax^ 

IT 



& 



{xo,yo,zo) 
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li{xo,yo,Zo + Az) 



oo 



E 



— (xo,yo,^o) 



s(a;o,|/o,^o) +E»=i {£1(2^0,2/0,^0)}] 



(xo,yo,^o) 



5x 



lim — — 

Ax^o Aa; 



{xq + Ax, 7/0, ^0) - h^x - Ax, 7/0, 2:0) 



(1) 



s (a;o,|/o,2:o) |^ 2Aa; 



s (xo, yo, ^0) + E»=i (0 (3^o> yo> ^0)}] 



dy 



lim — — 
Aj/->o Ay 



1 



(a^O) 1/0 + Ay, 2:0) - h^2' (a^o, yo - Ay, zq) 



(1) 



s (xo, yo, ^0) 



2Ay 



Ei^o ^ ||^ (^o,yo,2;o)| 



dn {x, y, z) 
dz 



lim — — 

Az->0 Az 



s(xo,yo,^o) + E»^i {0 (^o,yo,^o)} 

^3^^ (^^o, yo, ^0 + Az) - (xo, yo, ^0 - A2;) 



1 

s(xo,yo,2o) 



2Az 





"Ay» 1 
. 1 


Qyi+\ (a^o, yo, ^0) 






s(a;o,yo,^o) + Ej^i 


^¥{0(^0, 90,^0)}' 
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